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In order to popularize the so called Schwinger's method we reconsider the Feynman propagator of 
two non-relativistic systems: a charged particle in a uniform magnetic field and a charged harmonic 
oscillator in a uniform magnetic field. Instead of solving the Heisenberg equations for the position 
and the canonical momentum operators, R and P, we apply this method by solving the Heisenberg 
equations for the gauge invariant operators R and 7r = P — eA, the latter being the mechanical 
momentum operator. In our procedure we avoid fixing the gauge from the beginning and the result 
thus obtained shows explicitly the gauge dependence of the Feynman propagator. 

PACS numbers: 42.50.Dv 
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I. INTRODUCTION 

In a recent paper published in this journal [l|, three 
methods were used to compute the Feynman propaga- 
tors of a one-dimensional harmonic oscillator, with the 
purpose of allowing a student to compare the advan- 
tadges and disadvantadges of each method. The above 
mentioned methods were the following: the so called 
Schwinger's method (SM), the algebraic method and the 
path integral one. Though extremely powerful and el- 
egant, Schwinger's method is by far the less popular 
among them. The main purpose of the present paper is 
to popularize Schwinger's method providing the reader 
with two examples slightly more difficult than the har- 
monic oscillator case and whose solutions may serve as a 
preparation for attacking relativistic problems. In some 
sense, this paper is complementary to reference 

The method we shall be concerned with was in- 
troduced by Schwinger in 1951 2] in a paper about 
QED entitled "Gauge invariance and vacuum polariza- 
tion" . After introducing the proper time representa- 
tion for computing effetive actions in QED, Schwinger 
was faced with a kind of non-relativistic propagator in 
one extra dimension. The way he solved this problem 
is what we mean by Schwinger's method for comput- 
ing quantum propagators. For relativistic Green func- 
tions of charged particles under external electromagnetic 
fields, the main steps of this method are summarized 
in Itzykson and Zuber's textbook @ (apart, of course, 
from Schwinger's work Q). Since then, this method 
has been used mainly in relativistic quantum theory 
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However, as mentioned before, Schwinger's method is 
also well suited for computing non-relativistic propaga- 
tors, though it has rarely been used in this context. As 
far as we know, this method was used for the first time 
in non-relativistic quantum mechanics by Urrutia and 
Hernandez [171 ]. These authors used Schwinger's action 
principle to obtain the Feynman propagator for a damped 
harmonic oscillator with a time-dependent frequency un- 
der a time-dependent external force. Up to our knowl- 
edge, since then only a few papers have been written with 
this method, namely: in 1986, Urrutia and Manterola 
[l8T | used it in the problem of an anharmonic charged 
oscillator under a magnetic field; in the same year, Hor- 
ing, Cui, and Fiorenza [l|| applied Schwinger's method 
to obtain the Green function for crossed time-dependent 
electric and magnetic fields; the method was later applied 
in a rederivation of the Feynman propagator for a har- 
monic oscillator with a time-dependent frequency [2Cj| : 
a connection with the mid-point-rule for path integrals 
involving electromagnetic interactions was discussed in 
[2ll ] . Finally, pedagogical presentations of this method 
can be found in the recent publication [l[ as well as in 
Schwinger's original lecture notes recently published [Hj], 
which includes a discussion of the quantum action prin- 
ciple and a derivation of the method to calculate propa- 
gators with some examples. 

It is worth mentioning that this same method was inde- 
pendently developed by M. Goldberger and M. GellMann 
in the autumn of 1951 in connection with an unpublished 
paper about density matrix in statistical mechanics [23| . 

Our purpose in this paper is to provide the reader with 
two other examples of non-relativistic quantum propa- 
gators that can be computed in a straightforward way 
by Schwinger's method, namely: the propagator for a 
charged particle in a uniform magnetic field and this same 
problem with an additional harmonic oscillator potential. 
Though these problems have already been treated in the 
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context of the quantum action principle [181 ] , we decided 
to reconsider them for the following reasons: instead of 
solving the Heisenberg equations for the position and the 
canonical momentum operators, R and P, as is done in 
[HI], we apply Schwinger's method by solving the Heisen- 
berg equations for the gauge invariant operators R and 
7r = P — eA, the latter being the mechanical momen- 
tum operator. This is precisely the procedure followed 
by Schwinger in his seminal paper of gauge invariance 
and vacuum polarization This procedures has some 
nice properties. For instance, we are not obligued to 
choose a particular gauge at the beginning of calcula- 
tions. As a consequence, we end up with an expression 
for the propagator written in an arbitrary gauge. As a 
bonus, the transformation law for the propagator under 
gauge transformations can be readly obtained. 

In order to prepare the students to attack more com- 
plex problems, we solve the Heisenberg equations in ma- 
trix form, which is well suited for generalizations involv- 
ing Green functions of relativistic charged particles under 
the influence of electromagnetic fields (constant F^, a 
plane wave field or even combinations of both). For ped- 
agogical reasons, at the end of each calculation, we show 
how to extract the corresponding energy spectrum from 
the Feynman propagator. Although the way Schwniger's 
method must be applied to non-relativistic pro blems has 
already been explained in the literature (H Il8l. |22| , it is 
not of common knowledge so that we start this paper 
by summarizing its main steps. The paper is organized 
as follows: in the next section we review Schwinger's 
method, in section ITTT1 we present our examples and sec- 
tion [TV] is left for the final remarks. 



II. MAIN STEPS OF SCHWINGER'S METHOD 

For simplicity, consider a one-dimensional time- 
independent Hamiltonian Ti and the corresponding non- 
relativistic Feynman propagator defined as 



K {x, x'; t) = 9(t)(x\ exp 



-iHt 



(1) 



where 9(t) is the Heaviside step function and \x), \x') 
are the eingenkets of the position operator X (in the 
Schrodinger picture) with eingenvalues x and x' ', respec- 
tively. The extension for 3D systems is straightforward 
and will be done in the next section. For r > we have, 
from equation (jTJ) , that 



d 

ih—K(x. x ; t) 
Ot 



(x\Tiexp 



-iHt 



W). 



(2) 



Inserting the unity 11 = exp [— (i/T^Wr] exp [(i//l)7ir] in 
the r.h.s. of the above expression and using the well 
known relation between operators in the Heisenberg and 
Schrodinger pictures, we get the equation for the Feyn- 
propagator in the Heisenberg picture, 
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where \x,t) and |x',0) are the eingenvectors of opera- 
tors X(t) and A(0), respectively, with the correspond- 
ing eingenvalues x and x': X(t)\x,t) — x\x,t) and 
X(0)|a/,0) =x'\x',0), withK(x,x'\T) = (x, t\x' , 0). Be- 
sides, X(t) and P(r) satisfy the Heisenberg equations, 

m^(r) = [X(T) ) H} ; ift^(r) = [P(r),W]. (4) 
or ar 

Schwinger's method consists in the following steps: 

(i) we solve the Heisenberg equations for X(t) and P(r), 

and write the solution for P(0) only in terms of the 
operators X(t) and X(0); 

(ii) then, we substitute the results obtained in (i) into 

the expression for Ti(X(0), P(0)) in ^ and using 
the commutator [X(0),X(t)] we rewrite each term 
of Ti in a time ordered form with all operators X(t) 
to the left and all operators X(0) to the right; 

(iii) with such an ordered hamiltonian, equation can 
be readly cast into the form 

d 

ih—K(x 1 x'; t) = F(x, x'\ t)K(x, x';t), (5) 

OT 

with F(x,x';t) being an ordinary function defined 

as 

P , / ^_ {x,T\n ord (x(T),x(o))\ x ',o) 

[X,X ' >- {x,t\x',0) ■ W 

Integrating in r, the Feynman propagator takes the 
form 



K(x,x';t) = C(x,x')exp^~J F(x,x';T')dT ' j , (7) 

where C(x,x / ) is an integration constant indepen- 
dent of t and J means an indefinite integral; 

(iv) last step is concerned with the evaluation of 
C(x,x'). This is done after imposing the follow- 
ing conditions 

-ih-^(x,T\x\0) = (x,t\P(t)\x',0), (8) 
ih-^(x,r\x',0) = (x,t\P(0)\x',0), (9) 



as well as the initial condition 



lim K{x, x'\ t) = 5{x — x) 

T^0 + 



(10) 



ih—-K(x, x': t) 
Ot 



(x,t\H(X(0),P(0))\x',0), (3) 



Imposing conditions (jSJ) and ^ means to substitute in 
their left hand sides the expression for (x,t\x',0) given 
by 0, while in their right hand sides the operators P(r) 
and P(0), respectively, written in terms of the operators 
X(t) and X(0) with the appropriate time ordering. 
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III. EXAMPLES 

A. Charged particle in an uniform magnetic field 

As our first example, we consider the propagator of a 
non-relativistic particle with electric charge e and mass 
m, submitted to a constant and uniform magnetic field B. 
Even though this is a genuine three-dimensional problem, 
the extension of the results reviewed in the last section 
to this case is straightforward. Since there is no electric 
field present, the hamiltonian can be written as 



n = 



(P-eA) z 
2m 



7T 

2m ' 



(11) 



where P is the canonical momentum operator, A is the 
vector potential and ir = P — eA is the gauge invariant 
mechanical momentum operation. We choose the axis 
such that the magnetic field is given by B = Be 3 . Hence, 
the hamiltonian can be decomposed as 



H 



2m 



P 2 
2m 



P 2 
2m 



(12) 



with an obvious definition for 7i±. 

Since the motion along the OX 3 direction is free, 
the three-dimensional propagator K(x., x';r) can be 
written as a product of a two-dimensional propagator, 
K±(r, r'; r), related to the magnetic field and a one- 
dimensional free propagator, (x 3 , x' 3 ; r): 

K(x,x';r) = K x (v,T';T)Ki 0) (x 3 ,x' 3 ;T), (r > 0) 

(13) 

where r = x%ei + X2G.2 and (x 3 , a:!,; r) is the well 
known propagator of the free particle [24j | , 



Kf\x 3 ,x' 3 ]T 



2-Kihr 



exp 



im (x 3 — x'^y 



2h 



(14) 



In order to use Schwinger's method to compute the two- 
dimensional propagator K±(r, r';r) = (r,T|r',0), we 
start by writing the differential equation 

in^-(r,r\r',0) = (r,r|Wx(R±(0),irx(Q))|r , ,0> , (15) 

where Ri(r) = Xi(r)ei + X2(r)e2 and n±(r) — 
ni(T)e 1 +TT2(T)e2- In (fT5)) |r,r) and |r',0) are the eigen- 
vectors of position operators R(t) = X 1 (r)e 1 + X 2 (r)e 2 
and R(0) = X 1 (0)e 1 + X 2 (0)e 2 , respectively. More 
especifically, operators Xi(0), Xl(t), ^2(0) and X 2 (t) 
have the eigenvalues X -y^ X\ ^ X 2 and X2, respectively. 
In order to solve the Heisenberg equations for operators 
R^(t) and 7Tj_(t), we need the commutators 



= 2ihlli(T) , 

= 2iheBe ij3 ir j (T), (16) 



where 6ij 3 is the usual Levi-Civita symbol. Introducing 
the matrix notation 



Xi(t) 



7ri(r) 



R < T > - ( x~aA ) ■ n w = { I'M ) • < 17 > 

and using the previous commutators the Heisenberg 
equations of motion can be cast into the form 



rfR(r) _ n(r) 

dr = 
dn(r) 



dr 



m 

2wCTI(t) , 



(18) 
(19) 



where 2lo — eB/m is the cyclotron frequency and we 
defined the anti-diagonal matrix 



C = 



1 

-1 



Integrating equation (fTi?)) we find 

n(r) = e 2ujCT U(0) 



(20) 



(21) 



Substituting this solution in equation (|18[) and integrat- 
ing once more, we get 

R(r) - R(0) = ^i^ e " Cr n(0) , (22) 
mw 

where we used the following properties of C matrix: 
C 2 = -A; C~ 1 = -C = C T , e aC = cos (a) 11 + sin (a)C 
with C T being the transpose of C. Combining equations 
(|2"2"| and (|2"Tj) we can write 11(0) in terms of the operators 
R(t) and R(0) as 



n(0) = - 



sin (lot) 



R(t) - R(0) 



(23) 



In order to express fi± = (nf + it 2 )/2m in terms of R(r) 
and R(0), we use (|23p . In matrix notation, we have 



Hx = ^n T (o)n(o) 



mui 



2 sin (lot) 



R T (r)R(r) +R T (0)R(0) 



-R T (r)R(0) - R T (0)R(r) . (24) 



Last term on the r.h.s. of ([24]) is not ordered appropri- 
ately as required in the step (ii). The correct ordering 
may be obtained as follows: first, we write 

2 

R(0) T R(r) = R(r) T R(0) + ^[X(0), X t (r)] . (25) 



Using equation \22\ . the usual commutator 
[Xi(0), 7Tj(0)] = ihSijSL and the properties of matrix 
C it is easy to show that 



53^(0), *<(r)] 



2iftsin(a;r) cos(wr) 



mw 



(26) 
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so that hamiltonian Ti± with the appropriate time order- 
ing takes the form 



Ha 



2 sin 2 (ut) 
— ihco cot(wr) 



R 2 (r) + R 2 (0) - 2R T (r)R(0) 



(27) 



Substituting this hamiltonian into equation (|15[) and in- 
tegrating in t, we obtain 

, / < C(r,r') { imijj , . , 
(r,r|r',0) = , / { exp\ — — cot(wr)(r - r') 2 k (28) 



sin (ur 



2ft 



where C(r,r') is an integration constant to be deter- 
mined by conditions ((SJ , © and (fT0|) , which for the case 
of hand read 



,r|^(r)|r',0) 



(r,r|^-(0)|r',0) 



ih£--eMr)) (r.rr'.U) (20) 



iftA-eA^rOj (r,r|r',0> ,(30) 



lim (r,T|r',0) = c5 (2) (r - r'). 

T— >0+ 



(31) 



In order to compute the matrix element on the l.h.s. 
of (j29|) , we need to express II(r) in terms of R(r) and 
R(0). From equaitons ([2lTl and ([23]) . we have 

n(r) = _^ e ^ C f R(r) - R(0)) , (32) 
sm (ut) \ / 

which leads to the matrix element 

(r, r|7Tj(r)|r', 0} = mui[cot(uiT) {xj — x'A 

+ e jk3 (x k -x' k )){T,T\r',0) , (33) 

where we used the properties of matrix C and Einstein 
convention for repeated indices is summed. Analogously, 
the l.h.s. of equation (|3"0|) can be computed from (|23l . 

(r, T|7Tj(0)|r',0) = mui[cot(u>T) (Xj — x'j) 

- e jk3 (x k -x' k )]{r,T\r',0) . (34) 



Substituting equations (|3"3"j) and (031) into and 
(1301 . respectively, and using (f2"8")) . we have 



ia^- + e^(r)+ ^eF jk (x k -x' k )\c(v,v') = 0,(35) 
ift^T -eA,(r')+ ^fe^fc -4)] C(r,r 0=0,(36) 

where we defined F^-fc = 6^3 £?. 

Our strategy to solve the above system of differential 
equations is the following: we first equation (|35[) assum- 
ing in this equation variables r ' as constants. Then, we 
impose that the result thus obtained is a solution of equa- 
tion ([36]) . With this goal, we multiply both sides of (|35]) 
by dxj and sum over j, to obtain 



1 fdC 



le 
J 



Integration of the previous equation leads to 
C(r r ') = C(r ' r ') e { * r f r> lA ^^ )+ i F i* 

(38) 

where the line integral is assumed to be along curve T, to 
be specified in a moment. As we shall see, this line inte- 
gral does not depend on the curve T joining r ' and r, as 
expected, since the l.h.s. of ([37)) is an exact differencial. 

In order to determine the differential equation for 
C(r',r') we must substitue expression (|3"8"]) into equa- 
tion (|36p . Doing that and using carefully the fundamen- 
tal theorem of differential calculus, it is straightforward 
to show that 



— (r' r') =0 



(39) 



which means that C(r',r') is a constant, Co, indepen- 
dent of r '. Noting that 



[Bx(£-r% = -F jk (& -4) 
equation ([38]) can be written as 



(40) 



<7(r,r')=Co exp|| J *[A«) - ± B x (£ - r ')] -^j 

(41) 

Observe, now, that the integrand in the previous equa- 
tion has a vanishing curl, 



x 



A«)-±Bx(£-r') 



B = , 



which means that the line integral in (|42p is path in- 
dependent. Choosing, for convenience, the straightlinc 
from r ' to r, it can be readly shown that 



y r [Bx(|-r')M£ = 0, 



where r s ; means a straightline from r ' to r. With this 
simplification, the C(r',r) takes the form 



C(r,r')=C exp<j ^ r y A(£) • d$ 



(42) 



Substituting last equation into (|28|) and using the initial 
condition (fT0|) . we readly obtain Co = ^fn- Therefore 
the complete Feynman propagator for a charged particle 
under the influence of a constant and uniform magnetic 
field takes the form 



#(x,x';r) 



muj 



dxj . (37) exp 



2nih sin (lut 
irriLU 



2nihT 



2/; 
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where in the above equation we omitted the symbol r s ; 
but, of course, it is implicit that the line integral must be 
done along a straightline, and we brought back the free 
propagation along the OX3 direction. A few comments 
about the above result are in order. 

1. Firstly, we should emphasize that the line integral 
which appears in the first exponencial on the r.h.s. 
of (T4"3"l) must be evaluated along a straight line be- 
tween r' and r. If for some reason we want to choose 
another path, instead of integral J , A(£) ■ d£, we 
must evaluate / r r ,[A(£) - (1/2)B x (£ - r ')] • d£. 

2. Since we solved the Heisenberg equations for the 
gauge invariant operators Rj^ and tt±, our final 
result is written for a generic gauge. Note that 
the gauge-independent and gauge-dependent parts 
of the propagator are clearly separated. The gauge 
fixing corresponds to choose a particular expression 
for A(£). Besides, from (|43[1 we imediately obtain 
the transformation law for the propagator under a 
gauge transformation A — * A + VA, namely, 

K (r : r >■ T ) _ e f A(r) R ^ r /. r) e -fA(r ') _ 



ing this kind or regularization, we write 

/L/2 pL/2 L 2 mLO 
dxi I dx 2 K (r, r; -ih/3) « — 
-l/2 J -l/2 27rftsenh(fi/3w) 

L 2 eB 
2irh (e hPuJ - e - **") 

L 2 eB e -^ hf3wc 

_ ^ L eB p(n + l)fju c 
- 2^^H e 

n=0 

where we denoted by lu c = eB/2m the ciclotron 
frequency. Comparing this result with that of a 
partition function whose energy level E n has de- 
generacy g n , given by 

we imediately identify the so called Landau leves 
and the corresponding degeneracy per unit area, 

* -(„+!)*«., (n = 0,l,...). 



Although this transformation law was obtained in 
a particular case, it can be shown that it is quite 
general. 



3. It is interesting to show how the energy spectrum 
(Landau levels), with the corresponding degener- 
acy per unit area, can be extracted from prop- 
agator P3"| . With this purpose, we recall that 
the partition function can be obtained from the 
Feynman propagator by taking r = —ih(3, with 
(3 = 1/(KbT), and taking the spatial trace, 



dx\ I dxi K{r, r; —ih/3) . 



Substituting (143)) into last expression, we get 



Z{f3) 



dx\ I dx2 — 

-00 J —00 2nhseim(npuj) 



where we used the fact that sin(— iff) = —i sinh 8. 
Observe that the above result is divergent, since 
the area of the OX \X 2 plane is infinite. This is a 
consequence of the fact that each Landau level is 
infinitely degenerated, though the degeneracy per 
unit area is finite. In order to proceed, let us as- 
sume an area as big as we want, but finite. Adopt- 



B. Charged harmonic oscillator in a uniform 
magnetic held 

In this section we consider a particle with mass m and 
charge e in the presence of a constant and uniform mag- 
netic field B = Bea and submitted to a 2-dimensional 
isotropic harmonic oscillator potential in the OX\X 2 
plane, with natural frequency luq. Using the same no- 
tation as before, we can write the hamiltonian of the 
system in the form 



p2 

2m 



where 



n± = 



2m 



—muj 
2 



Xn 



(44) 



(45) 



As before, the Feynman propagator for this problem 
takes the form AT(x,x';r) = K±(r, r'; t)K^ (23, x' 3 ; r), 
with K$(x3, x' 3 ; r) given by equation (fT4|) . The propa- 
gator in the OA'iA^-plane satisfies the differential equa- 
tion (fT5|) and will be determined by the same used in the 
previous example. 

Using hamiltonian (|4l))) and the usual commutation re- 
lations the Heisenberg equations are given by 

dR(r) n(r) 



dr 
cgl(r) 

dr 



2uCU(t) - to^R(t) , 



(46) 
(47) 
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where we have used the matrix notation introduced in 
(fTTj) and (HOI)- Equation gB]) is the same as HU), but 
equation (|47| contains an extra term when compared to 
(Unj). In order to decouple equations ([4*6")) and (|4"T)) . we 
differentiate ([46]) with respect to r and then use (|47|) . 
This procedure leads to the following uncoupled equation 

d 2 R(r) dRM 2 , , 

, o - 2wC— + w 2 R(r) = (48) 

CIT CIT 

After solving this equation, R(r) and II(r) are con- 
strained to satisfy equations (f4"6"]) and (|4"7]) , respectively. 
A straightforward algebra yields the solution 



R(r) = MTR(0) + NII(0) 
n(r) = M+II(0) - m 2 w 2 NR(0) 



where we defined the matrices 
sin (fir) , „ip 



mfi 



cos (fir) 11 ± -sin(fir)C 



(49) 
(50) 



(51) 
(52) 



and frequency fi = y/uj 2 + uj 2 . Using (|49|) and (|50|) . we 
write 11(0) and II(r) in terms of R(r) and R(0), 



n(0) = N _1 R(r) - N _1 M _ R(0) , 



(53) 



n(r) = 



"N^Rfr)- 



N~ M"+m WqN R(0). (54) 



Now, we must order appropriately the hamiltonian 
operator H±_ = II T (0)n(0)/(2m) + mw 2 R T (0)R(0)/2, 
which, with the aid of equation (|53|) . can be written as 



Ha 



R t (t)(N- 1 ) t - R T (0)(M-) T (N- 1 ) T 
N -1 R(r) - N _1 M _ R(0) + mw 2 R T (0)R(0) 

R T (r) -R T (0)(M-) 3 



2 sin 2 (fir) 



R(r) - M-R(O) + mw 2 R T (0)R(0) 



2 sin 2 (fir) 



R T (r)R(r) - R T (r) 



-R(0) 



-R T (0)(M") T R T (r) + R T (0)(MT) T MTR(0) 
+mw 2 R T (0)R(0) 



2 sin 2 (fir) 



R 2 (r) -R T (r)M~R(0) 
R T (0)(MT) T R(r) +R 2 (0) 



(55) 



where superscript T means transpose and we have used 
the properties of the matrices N and M~ given by (|5"Tj) 
and ([52]) . In order to get the right time ordering, observe 
first that 



where 
(M-RCOjV.XiCr) 



zfi.Tr 



— - sin (2fir) . 
mfi 



Using the last two equations into f|55|) we rewrite the 
hamiltonian in the desired ordered form, namely, 

mfi 2 



2 sin 2 (fir) 



BT(t) + R 2 (0) - 2R 1 (t)M-R(O) 



sin (2fir 

mS2 



For future convenience, let us define 



(56) 



U(t) = cos (lut) cos (fir) + — sin (wt) sin (fir) , (57) 
V(t) = sin (ur) cos (fir) — — cos (wt) sin (fir) (58) 



and write matrix M , defined in (|52p . in the form 

M~ = J7(r)l + V(r)C. (59) 
Substituting ([59]) in |56|) we have 
mfi 2 



n± = 



2 sin 2 (fir) 



R 2 (r) + R 2 (0) - 2£/(r)R J (r)R(O) 



/'ft 



2U(r)R J (r)CR(O) sin(2fir) 

mfi 



(60) 



The next step is to compute the classical function 
F(r, r'; r). Using the following identities 



fit/(r) 
sin 2 (fir) 

ny(r) 

sin 2 (fir) 



d 
"dr 
d 

~ch 



cos (wt) 



sin (fir) 
sin (wr) 



sin (fir) 



(61) 
(62) 



into (JHUJ) , we write F(r, r';r) in the convenient form 
F(r,r';r) 

" sin (wt) 



mfi 2 , 2 ,2 2 ^ . drcos(wr) 

(r 2 + r' )csc(fir) 2 + mfir • r' ) — f 

2 drLsin(fir) 

cos (fir) 



+ mfir • Cr' 



fir 



sin (fir) 



ifi.fi 



sin (fir) 



(63) 



Inserting this result into the differential equation 

jfiA(r,r|r',0) = F(r, r'; r)(r, r|r', 0) , 

and integrating in r, we obtain 

/ i / n\ C '( r ' r ') f 
(r, r|r ,0) = — — 77r-r ex P 



sin (fir) I 2fi 
,cos (wt) 



(r 2 + r' ) cot (fir) 



2 r r 



sin (fir) 



r-Cr' 



sin (fir) 



(64) 



R T (0)(M~) T R(r) = R T (r)M"R(0)- 



where C(r, r') is an arbitrary integration constantto be 
determined by conditions $2§j), §U\) and Using (JSl) 
we can calculate the l.h.s. of condition ([29"]) . 

TTifi f 

(r,T|7rj(T)|r',0) = sin Fq3 { cos ( r?T ) a: 3 -cos(wr)^- 

-sin(fir)a; fc -sm{uT)x' k e jfe3 |(r, r|r', 0), (65) 



7 



and using (153)) we get the l.h.s. of condition ((30)) , 

(r,r[7r,-(0)|r',0) = - — cos (u;r)x 7 - - cos (fir) x', 
sin (Ut) I J 

+ ^ sin (f2r)xfc - sin (wr)iTfc 6jfe3}(r,r|r',0}. (66) 
With the help of the simple identities 



where the path of integration T will be specified in a 
moment. Inserting expression (|7ip into the second differ- 
ential equation (|70|) . we get 



C(r',r')=C ., 



A (r2 + / ) = 2x . . _^ (r2 + / ) = 2 ^ 



<9;C 



- — r • r = x ■ 



r • r = Xj 



d d 
—r-Cr=e jk3 x k ; —r-Cr' = -e jk3 x k . 



and also using equation (|64p . we are able to compute the 
right hand sides of conditions (|29f and (|30|) . which are 
given, respectively, by 



ih dC(r,r') n cos (fir) 
1 ml 2 



C(r,r') as 3 - 



Cj — mfi 



cos (cut) , 



sin (fir) ~ J sin (fir) J 

mfi^&4e jfe3 4 " e^(r)}(r,r|r',0) (67) 



sin (fir) 



where Co is a constant independent of r ', so that equa- 
tion (|71[) can be cast, after some convenient rearrange- 
ments, into the form 



C(r,r') =C exp 



r. 



A(£) - x C 



(72) 

Note that the integrand has a vanishing curl so that we 
can choose the path of integration T at our will. Choos- 
ing, as before, the straight line between r' and r, it can 
be shown that 



A(£)-fx£ 



= £ A(0'^+^r-Cr', (73) 



and 



ih dC(r,r') cos(f2T) , cos(wr) 
C(r,r') dx] sin (Qt) ^ + sin (fir) Xj 

sin (Or) 



0k3 x k -eA 3 {r')Yr,T\Y\Q). (68) 



Equating 455} and 457} , and also 455} and 155}), we get 
the system of differential equations for C(r, r') 



in — + e 

axj 



C(r,r')=0, (69) 



A 3 (v) + ^fx k 

A 3 {v') + ^x' k C(r,r') = 0. (70) 



Proceeding as in the previous example, we first integrate 
(|69[) . With this goal, we multiply it by dxj, sum in j and 
integrate it to obtain 



(71) 



where, for simplicity of notation, we omitted the symbol 
T s ; indicating that the line integral must be done along 
a straight line. From equations ([71"} . ([72} e ([73} . we get 



C(r,r')=C expi| / A(£) • [ exp | • Cr'^ . 



which substituted back into equation (|64|) yields 



(74) 



eX K2^W{ (r2+r ' 2)COS(fiT) 
-2r • r' cos (ut) - 2 sin (lot) - ^ sin VLt rCr'jj (75) 



The initial condition implies Co = mf2/ '(2mh). Hence, 
the desired Feynman propagator is finally given by 



if(x,x';r) = K ± (r,r';T)Ki \x 3 ,x' 3 ;T 



ex P \ -t~ I M£)-d£ }e*P 



2 7r i h sin (Or) V 27riftr 



2 cos (wr)r ■ r' — 2 sin (wt) — — sin (fir) 



2hsm(Qr) 



|cos (Or)(r 2 + r' 2 ) 



4} 



exp 



im (x 3 — S3) 2 
2fr r 



(76) 
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where we brought back the free part of the propagator 
corresponding to the movement along the OX 3 direction. 
Of course, for w = we reobtain the propagator found 
in our first example and for B = we reobtain the prop- 
agator for a bidimensional oscillator in the OX 1X2 plane 
multiplied by a free propagator in the OX3 direction, as 
can be easily checked. 

Regarding the gauge dependence of the propagator, 



the same comments done before are still valid here, 
namely, the above expression is written for a generic 
gauge, the transformation law for the propagator under a 
gauge transformation is the same as before, etc. We fin- 
ish this section, extracting from the previous propagator, 
the corresponding energy spectrum. With this purpose, 
we first compute the trace of the propagator, 



O />OC 

dxi dx2Kj_(xi,xi,X2,xa;T) 

x: J —00 



in 



2nih sin(Or) 
1 



dxi dx2 exp 



2h sin(nr) 



jfcos(Qr) — cos(wr)^ (x\ + x^) 



2[cos(S1t) — cos(wt)] 



(77) 



r 



where we used the well known result for the Fresnel in- 
tegral. Using now the identity 

cos(rJr)-cos(wr) = -2 sin[(f2+£j)r/2] sin[(fi-o;)r/2)] , 
we get for the corresponding energy Green function 

/>oc />oo />oo 

Q(E)=—i dr e^ ST / dx\ I dx 2 K±(xi, Xi, x 2 , x 2 ; r) 

JO J — 00 J — 00 



IV. FINAL REMARKS 



dr- 



en 



Et 



sen(^fr) sen(^r) 



poo ( 00 



where is tacitly assumed that E — > E — is and we also 
used that (with the assumption v ^ v — ie) 



-i(ra+i)i. 



sen(^r) 



n— 



Changing the order of integration and summations, and 
integrating in r, we finally obtain 



l,n=0 



E — E n i 



(78) 



where the poles of Q(E), which give the desired energy 
levels, are identified as 

E nl = (l + n + l)m + (l-n)tko, (l,n = 0, 1, ...) . (79) 

The Landau levels can be reobtained from the previous 
result by simply taking the limit wq — > 0: 

E nl — » (21 + l)tkj = (l+ i)fiw c , (80) 

with I = 0,1,... and lj c — eB/m, in agreement to the 
result we had already obtained before. 



In this paper we reconsidered, in the context of 
Schwinger's method, the Feynman propagators of two 
well known problems, namely, a charged particle under 
the influence of a constant and uniform magnetic field 
(Landau problem) and the same problem in which we 
added a bidimensional harmonic oscillator potential. Al- 
' though these problems had already been treated from 
the point of view of Schwinger's action principle, the 
novelty of our work relies on the fact that we solved 
the Heisenberg equations for gauge invariant operators. 
This procedure has some nice properties, as for instance: 
(i) the Feynman propagator is obtained in a generic 
gauge; (ii) the gauge-dependent and gauge-independent 
parts of the propagator appear clearly separated and 
(Hi) the transformation law for the propagator under 
gauge transformation can be readly obtained. Besides, 
we adopted a matrix notation which can be straightfor- 
wardly generalized to cases of relativistic charged parti- 
cles in the presence of constant electromagnetic fields and 
a plane wave electromagnetic field, treated by Schwinger 
[2j. For completeness, we showed explicitly how one can 
obtain the energy spectrum directly from que Feynman 
propagator. In the Landau problem, we obtained the (in- 
finitely degenerated) Landau levels with the correspond- 
ing degeneracy per unit area. For the case where we 
included the bidimensional harmonic potential, we ob- 
tained the energy spectrum after identifying the poles of 
the corresponding energy Green function. We hope that 
this pedagogical paper may be useful for undergraduate 
as well as graduate students and that these two simple 
examples may enlarge the (up to now) small list of non- 
relativistic problems that have been treated by such a 
powerful and elegant method. 



9 



Acknowledgments 

F.A. Baronc, H. Boschi-Filho and C. Farina would like 
to thank Professor Marvin Goldberger for a private com- 



munication and for kindly sending his lecture notes on 
quantum mechanics where this method was explicitly 
used. We would like to thank CNPq and Fapesp (brazil- 
ian agencies) for partial financial support. 



[1] F.A. Barone, H. Boschi-Filho and C. Farina, Three meth- 
ods for calculating the Feynrnan propagator, Am. J. Phys. 
71, 483-491 (2003). 

[2] J. Schwinger, On Gauge Invariance And Vacuum Polar- 
ization, Phys. Rev. 82, 664 (1951). 

[3] Claude Itzykson and Jean-Bernard Zuber, Quantum 
Field Theory, (McGraw-Hill Inc., NY, 1980), pg 100. 

[4] E.S Fradkin, D.M Gitman and S.M. Shvartsman, Quan- 
tum Eletrodinamics with Unstable Vacuum (Springer, 
Berlim, 1991). 

[5] V. V. Dodonov, I. A. Malkin, and V. I. Manko, "In- 
variants and Green's functions of a relativistic charged 
particle in electromagnetic fields," Lett. Nuovo Cimento 
Soc. Ital. Fis. 14, 241-244 (1975). 

[6] V. V. Dodonov, I. A. Malkin, and V. I. Manko, "Green's 
functions for relativistic particles in non-uniform external 
fields," J. Phys. A 9, 1791- 1796 (1976). 

[7] J. D. Likkcn, J. Sonnenschein, and N. Weiss, "The theory 
of anyonic superconductivity: A review," Int. J. Mod. 
Phys. A 6, 5155-5214 (1991). 

[8] A. Ferrando and V. Vento, "Hadron correlators and the 
structure of the quark propagator," Z. Phys. C63, 485 
(1994). 

[9] H. Boschi-Filho, C. Farina and A. N. Vaidya, 
"Schwinger's method for the electron propagator in a 
plane wave field revisited," Phys. Lett. A 215, 109-112 
(1996). 

[10] S. P. Gavrilov, D. M. Gitman and A. E. Goncalves, "QED 
in external field with space-time uniform invariants: Ex- 
act solutions," J. Math. Phys. 39, 3547 (1998). 

[11] D. G. C. McKeon, I. Sachs and I. A. Shovkovy, "SU(2) 
Yang-Mills theory with extended supersymmetry in a 
background magnetic field," Phys. Rev. D59, 105010 
(1999). 

[12] T. K. Chyi, C. W. Hwang, W. F. Kao, G. L. Lin, 
K. W. Ng and J. J. Tseng, "The weak-field expansion for 
processes in a homogeneous background magnetic field," 
Phys. Rev. D 62, 105014 (2000). 

[13] N. C. Tsamis and R. P. Woodard, "Schwinger's propaga- 
tor is only a Green's function," Class. Quant. Grav. 18, 



83 (2001). 

[14] M. Chaichian, W. F. Chen and R. Gonzalez Felipe, 
"Radiatively induced Lorentz and CPT violation in 
Schwinger constant field approximation," Phys. Lett. B 
503, 215 (2001). 

[15] J. M. Chung and B. K. Chung, "Induced Lorentz- 
and CPT-violating Chern-Simons term in QED: Fock- 
Schwingcr proper time method," Phys. Rev. D 63, 
105015 (2001). 

[16] H. Boschi-Filho, C. Farina e A.N. Vaidya, Schwinger's 
method for the electron propagator in a plane wave field 
revisited, Phys. Let. A215 (1996) 109-112. 

[17] Luis F. Urrutia and Eduardo Hernandez, Calculation of 
the Propagator for a Time- Dependent Damped, Forced 
Harmonic Oscillator Using the Schwinger Action Princi- 
ple, Int. J. Theor. Phys. 23, 1105-1127 (1984) 

[18] L.F. Urrutia and C. Manterola, Propagator for the 
Anisotropic Three-Dimensional Charged Harmonic Os- 
cillator in a Constant Magnetic Field Using the 
Schwinger Action Principle, Int. J. Theor. Phys. 25, 75- 
88 (1986). 

[19] N. J. M. Horing, H. L. Cui, and G. Fiorenza, Nonrel- 
ativistic Schrodinger Greens function for crossed time- 
dependent electric and magnetic fields, Phys. Rev. A34, 
612615 (1986). 

[20] C. Farina and Antonio Segui-Santonja, Schwinger's 
method for a harmonic oscillator with a time- dependent 
frequency, Phys. Lett. A184, 23-28 (1993). 

[21] S.J. Rabello and C. Farina, Gauge invariance and the 
path integral, Phys. Rev. A51, 2614-2615 (1995). 

[22] J. Schwinger, Quantum Mechanics: Symbolism of Atomic 
Measurements, edited by B.G. Englert (Springer, 2001). 

[23] Private communication with Professor M. Goldberger. 
We thank him for kindly sending us a copy of his notes 
on quantum mechanics given at Princeton for more than 
ten years. 

[24] R.P. Feynman e A.R. Hibbs, Quantum Mechanics and 
Path Integrals (McGraw-Hill, New York, 1965). 



